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Abstract 

In this paper, we classify the perfect lattices in dimension 8. There are 10916 
of them. Our classification heavily relies on exploiting symmetry in polyhedral 
computations. Here we describe algorithms making the classification possible. 

1 Introduction 

A lattice L is a rank d subgroup of M*^. It is of the form L = TLv\ + ■ ■ ■ + Zw^ for 
linearly independent Ui, . . . , which are called a lattice basis of L. The determinant 
detL = I det(fi, . . . ,Vd)\ is independent of the chosen basis. By i?'' we denote the 
Euclidean unit ball. 

The packing radius A(L) of L is defined as half of the minimum distance between 
distinct lattice points. We call L + X{L)B'^ the lattice packing defined by L. The pack- 
ing density of the lattice packing defined by L is A(L)'^ vol(i?°')/ det L, the proportion 
of space covered by balls. The lattice packing problem consists in finding the lattices 
having highest packing density. 

A quadratic form q{x) over M'^ is a function defined as q{x) = ^xQx with Q a 
real [d x d) symmetric matrix. The quadratic form q is called positive definite if 
the corresponding matrix Q is positive definite. We denote by S'^ (resp. 5>o) 
set of symmetric matrices (resp. positive definite matrices). In this paper, we will 
identify forms with their symmetric matrices. Two matrices A,B & iS^q are called 
arithmetically equivalent if there exists a P G GLrf(Z) such that B = ^PAP. Denote 
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by Aut(A) the arithmetic automorphism group of A G 5^, i.e., the group of all 
P e GLrf(Z) such that A = 'PAP. 

If one takes a basis V = {vi, . . . ,Vd) of a lattice L, then L can be reconstructed, 
up to orthogonal transformations, from the Gram matrix Gy = {{vi^Vj)) G S'^q with 
( , ) being the standard scalar product on Mf'. Another basis V of L corresponds to 
another Gram matrix Gy'-, which is arithmetically equivalent to Gy- 

In other words, isometry classes of lattices are in one-to-one correspondence with 
arithmetical equivalence classes of positive definite quadratic forms. For computa- 
tions, it is preferable to use the language of positive definite quadratic forms, which 
we will do in the remaining of the paper. For more details on the correspondence 
between lattice properties and quadratic forms, see |Ma03l Chapter 1]. 

The arithmetical minimum of a positive definite form A is \{A) = min^,g2d_|o} *vAv. 
The Hermite invariant of A is defined as 7(A) = (j^^^^iyd- The packing density of 

the lattice corresponding to A is a/7(A)'^ ™2^ \ with vol(i?'^) being the volume of the 
unit ball B'^. Hence, solving the lattice packing problem is equivalent to maximizing 
7 over SyQ. Denote by 7^ = max^g^d^ 7(^) the Hermite constant in dimension d. 
The set of shortest vectors of A is 

Mm{A) = {veZ'^ I 'vAv = X{A)} . 

The form A is called perfect if the equations 

%Bv = X{A) for every v G Min(y4) 

have the unique solution B = A among B E S"^. Perfect forms were introduced in 
|KoZo73j . and studied by many authors, for example, in |KoZo77j . |Vo08j . |CoS188j . 
[Co5T] and |Ry70l . See also [Z^ and [Mi03] . 

A form A is called eutactic if there exist A^, > such that 

A-^ = J2 • 

ueMin(yl) 

A form A is called extreme if the Hermite invariant 7 attains a local maximum at A. 
Voronoi proved (see [V0O8I § 17] and |Ma03t Theorem 3.4.6]) that a form is extreme 
if and only if it is perfect and eutactic. 

Given a perfect form A, its perfect domain is the polyhedral cone 

Dom(A) = 1 ^ X^v'v I A„ > I . 

v&Mm{A) ^ 

Perfect domains form a face-to-face tessellation containing the cone iS^q, i.e., every 
Q G belongs to at least one perfect domain and the intersection of two perfect 
domains Domi, Dom2 is a face of both. If two perfect forms A, B are arithmetically 
equivalent, i.e., A = 'PBP for some P G GLrf(Z), then Dom(A) = P Dom(i?) *P. 
Voronoi proved (see |Vo08l § 7] and |Ma03l Theorem 7.4.4]) that for every fixed 
dimension d, the number of arithmetical inequivalent perfect forms is finite. Note, 
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Table 1: Perfect forms up to dimension 8. In fact, the lattices realizing the maximum 
of 7 are root lattices (see |CoS198l Chapter 6] or |Ma03i Chapter 4]) 

since perfect domains form a tessellation of iSJ^q, their interest is larger than only 
sphere packing theory (see, for example, |AMRT75] . |Sh05j and |So05] ). 

Since the group GLrf(Z) acts on the tessellation with perfect domains, this provides 
an algorithm for classifying all perfect forms in a fixed dimension. In particular one 
can solve the lattice packing problem using this algorithm. 

Perfect forms have been classified up to dimension 7: In [KoZo77j Korkine and 
Zolotarev obtained the classification of perfect forms up to dimension 5 without using 
Voronoi's algorithm. The classification of perfect forms in dimension 6, 7 was done 
by Barnes (see |Ba57j ) and Jaquet (see |Ja93] ) using Voronoi's algorithm. 

Theorem 1.1 There are 10916 perfect forms in dimension 8. 

We prove the above theorem by implementing Voronoi's algorithm (see |La71] ). This 
enumeration problem in dimension 8 was considered by Martinet and his school: in 
fact after the work of Lai'hem |La92j . Baril |Ba96j . Napias [Na96j and Batut and 
Martinet |BaMa05"] a list of 10916 perfect forms was known and our contribution 
consists in proving that this list is complete. One key step of the enumeration is to 
prove: 

Theorem 1.2 The polyhedral cone Dom(QE8) has 25075566937584 facets in 83092 
orbits. 

A direct consequence of the enumeration is: 

Theorem 1.3 (Conjecture 6.6.7 of lMa03\ l) Every perfect 8-dimensional lattice has 
a basis of minimal vectors. 

Using the face structure of the perfect domains we obtain: 

Theorem 1.4 The set of possible kissing numbers \ Min(y4)| for A G S^q is 2{1 . . . 58, 
60, 63... 71, 75, 120}. 

Using Theorem 11.11 Riener |Ri06] classified all extreme lattices in dimension 8: 
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Theorem 1.5 There are 2408 extreme lattices in dimension 8. 

Note that Mordell's inequality 7d+i < if^'^'^^ (for a proof see [Z^ 
Section 2.3] and |Ma03t Section 2.3]) and the knowledge of 77 and Eg yields the value 
of 78- The value of 73 was first computed by Blichfeldt in |B135j . The first proof 
that Eg is the unique form realizing 73 is in |Ve80j . Another proof based on analyzing 
the equality case in Mordell's inequality can be found in |Ma03t Section 6.6]. A 
completely different, computer assisted proof is in |CoKu04j . Our classification gives 
a fourth proof for the fact that Eg is the unique lattice attaining 73. 

Basic algorithms necessary for implementing Voronoi's algorithm are explained in 
Section [2l in Section [3] some dual description algorithms used are explained. In Sec- 
tion H] we explain specific implementation details and in Section [5] some consequences 
of our computation. 

2 Basic algorithms 

2.1 Notions from polyhedral combinatorics 

By the Farkas-Minkowski-Weyl Theorem (see e.g. |Sch86[ Corollary 7.1a]) a convex 
polyhedral cone C C M"^ is defined either by a finite set of generators {vi, . . . , vn} ^ 
or by a finite set of linear functionals {/i, . . . , /m} ^ (M™)*: 



C is called full- dimensional if the only vector space containing it is M™". C is called 
pointed if no linear subspace of positive dimension is contained in it. 

Let C be a full-dimensional pointed convex polyhedral cone in M™'. Given / G 
(M™)*, the inequality f{x) > is said to be valid for C if it holds for all x E C. A 
face of C is a pointed polyhedral cone {x E C \ f{x) = 0}, where f{x) > is a valid 
inequality. 

A face of dimension 1 is called an extreme ray of C; a face of dimension m — 1 is 
called a facet of C. The set of faces of C forms a partially ordered set under inclusion. 
We write F <iG if F (Z G and dim F = dim G — 1. Two extreme rays of C are said to 
be adjacent if they generate a two-dimensional face of C. Two facets of C are said to 
be adjacent if their intersection has dimension m — 2. Any (m — 2)-dimensional face 
of C is called a ridge and it is the intersection of exactly two facets of C. 

Every minimal set of generators {vi, . . . ,vn'} defining C has the property 



Every minimal set of linear functionals {/i, . . . , /m'} defining C has the property that 
{Fi, . . . , Fm'} with Fi = {x E C \ fi{x) = 0} is the set of facets of C. The problem of 
transforming a minimal set of generators into a minimal set of linear functionals (or 
vice versa) is called the dual description problem. 




i=l 



{M+fi, . . . , M+f AT/} = {e I e extreme ray of C}. 
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2.2 Voronoi's algorithm for classifying perfect forms 

In this section we describe Voronoi's algorithm. It computes a complete representative 
system of arithmetically inequivalent perfect forms: 

Input: Dimension d. 

Output: Set M. of all inequivalent d-dimensional perfect forms. 
7W ^ 0. 

while there is & Q eT do 
M^MU{Q}. 
T^T\{Q}. 
T *^ facets of Dom(Q). 
for F G do 

find perfect form Q' with F = Dom{Q) fl Dom{Q'). 
if Q' is not equivalent to a form in U T then 

T ^TU{Q'}. 
end if 
end for 
end while 

For the quadratic form Qa^ in the above algorithm, we may use (5a„ = {Qi,j)i<i,j<n 
with qi^i = 2, qi^i-i = qi-i/i = —1 and qij = otherwise (see [CoS198[ Section 6.1] or 
iMiQHl' Section 4.2]). 

By Voronoi's finiteness theorem, we know that the above program will eventually 
finish. 

The dual description part, which is used in computing the facets of Dom{Q), is 
computationally the most demanding part. The special methods used for that purpose 
are explained in Section 3. The computation of the adjacent domain is explained in 
Section 12.31 and the test of isometry of lattices in Section 12. 4[ 

2.3 Adjacent domain and shortest vector problems 

In this section we describe the subalgorithm which computes the adjacent perfect 
domain. Given a positive definite form A we need to solve the shortest vector problem, 
i.e., compute its arithmetical minimum X{A) and the set of vectors Min(A) that 
realize it. The Fincke-Pohst algorithm (cf. [Co93l Algorithm 2.7.7]), which has many 
implementations (sv |Va99j . in GAP, in MAGMA, in PARI, etc.), does this. In 
particular, solving a shortest vector problem in dimension 8 is a routine task and 
takes only a fraction of a second. 

If F denotes a facet of Dom(y4) with A a perfect form, then we define 

MmpiA) = {v e Mm{A) \ v^v e F} . 

We then have the following algorithm: 
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Input: Perfect form A and facet F of Dom(A). 
Output: Perfect form A' with Dom(A') n Dom(y4) = F. 

w ^ an element of Minp{A). 

U ^ solution in S*^ of ^vUv = for v G Min^(y4) and ^vUv = 1 for a 

V e Mm{A) - MmpiA). 
A ^ 1. 

while Min(A + \U) C Min(A) do 

A ^ 2A. 
end while 

while there is a f o G Z" — {0} such that ^Vo{A + XU)vo < *w(a4 + XU)w do 

A ^ solution of ^vo{A + XU)vo = ^w{A + XU)w. 
end while 
A' ^A + XU. 

How do we test if there is a fo G — {0} such that ^vqA'vq < ^wA'wl If A' is positive 
definite this is done by solving a shortest vector problem. If A' is positive semidefinite 
but not positive definite we can find a vector f o G Z" — {0} in the kernel of A' , which 
will satisfy ^VqA'vq = 0. If A' is not definite, we can find a vector f o G Z" — {0} such 
that '^VqA'vq < by taking rational approximations of an eigenvector corresponding 
to a negative eigenvalue. 

Note that any P G Aut(74) defines an action on Min(y4) by x i— > Px. This makes 
the group Aut(y4) act on Min(74), on the facets of Dom(y4) and on the perfect domains 
adjacent to Dom(A). 

2.4 Isometry tests 

In order to implement Voronoi's algorithm, we need to be able to decide if two perfect 
forms A and B are arithmetically equivalent or not. We also need for the dual 
description algorithm explained later a way to compute Aut(y4). 

If A and B are perfect forms, then B = ^PAP if and only if PMin(y4) = Min(i?). 
One implication is immediate; the other follows from the definition of perfect forms. 
Hence a possible algorithm is to compute Min(y4) and Min(i?) and search for a P G 
GL„(Z) such that PMin(y4) = Min(P); the search space is finite since Min(A) and 
Min(P) are finite and form a generating system of R". 

For a non-perfect form the set of minimal vectors does not characterize it. How- 
ever, a form A E S'^ is uniquely characterized by the set of all vectors f G Z*^ 
with norm ^vAv less or equal to the maximum diagonal coefficient ma.xf^-^ Ai^i. The 
program Isometry (see |P10pSc98| , |PlSo97] ). and its companion Aut_Grp for com- 
puting the automorphism group of a form, work by generating these possibly very 
large vector sets. Hence in order to speed up computations with these programs, it 
is desirable to compute with forms with small, maximum diagonal coefficients. This 
can be achieved by storing only a Minkowski reduced form for each equivalence class. 
Such can be obtained for example with the program Mink_red (see |P10pSc9"8| ). 
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3 Dual description methods 



General purpose programs like cdd jFii95 ]. Irs |Av93j . pd |Ma97j and porta |ChLo97] 
allow one to compute the dual description of a polyhedral cone given by its facets 
(linear inequalities) or by its extreme rays (generators). Since the programs are 
implementations of quite different methods, their efficiency may vary tremendously on 
a particular cone. The perfect form Qes (see |CoS198l Chapter 6]) has | Min((5Eg)| = 
240, therefore the 36-dimensional cone Dom(QE8) has 120 extreme rays. All general 
purpose programs take too much time in computing the facets of this cone. 

Usually one is interested only in a list of representatives of orbits of facets and 
not in the full list of facets. This leads us naturally to the Adjacency Decomposition 
Method, which exploits the symmetry of a polyhedral cone. The programs mentioned 
above are still used, but as a subroutine. In the following discussion, we will assume 
that C is a full dimensional, pointed polyhedral cone in M"^ generated by the extreme 
rays (ei)i<i<Ar and we want to compute its facets. 

3.1 Adjacency Decomposition Method 

Input: Extreme rays of a polyhedral cone C and a group G acting on C. 
Output: Complete set M. of inequivalent facets of C under G. 

T ^ {F} with F a facet of C. 

while there is a F G T do 
M ^MU{F}. 
T ^T\{F}. 
T ^ facets of F. 
ior H eJ" do 

find facet F' of C with H = FnF'. 

if F' is not equivalent under G to a facet in U T then 

T ^TU{F'}. 
end if 
end for 
end while 

Note that Voronoi's algorithm is very similar to the Adjacency Decomposition Method. 
Both fit into the framework of graph traversal algorithms. In the above algorithm, an 
initial facet can be found by solving a linear program. 

The algorithm relies on the ability to test if two facets are equivalent under the 
symmetry group. The possible strategies and the one that we used are explained in 
Section HI The Adjacency Decomposition Method framework is a reasonably natu- 
ral algorithm for computing with symmetry. Hence it was discovered several times, 
for example, in |Ja93j as "algorithm de I'explorateur" , in |ChRe96j as "adjacency 
decomposition method" and in [DFPSOl] as "subpolytope algorithm". 

We explain here the gift-wrapping step (see |Swa85j ) to compute an adjacent facet. 
Given the list of extreme rays (ej)i<j<Ar, a facet F <\ C is encoded by an index set 



7 



Sp C {I, . . . , N} such that F is generated by {ei)i^sF- Given a ridge H, we need an 
algorithm for computing the uniquely determined facet F' of C such that FdF' = H. 
The ridge H is encoded by a set Sh C Sp such that {ei)i^sH generate H. The defining 
inequalities / G (R™)* of the facet F' should satisfy /(e,) = for all i G Sh- The 
vector space of such functions has dimension 2. Let us select a basis {/i,/2} of it. 
If / = ai/i + 0^2/2 is the defining inequality of F or F', /(cj) > for all i with 
1 < i < N . This translates into a set of linear inequalities on ai,a2 defining a 2- 
dimensional pointed polyhedral cone. One easily finds its two generators {a\, a2)i<j<2- 
The corresponding inequalities fi{x) = a\fi{x) + 02/2(2;) > on C define the two 
adjacent facets F and F' of C. 

The Adjacency Decomposition Method can find the dual description of very sym- 
metric polyhedral cones, when other methods fail. But this algorithm uses dual 
description, albeit in one dimension lower and again this computation might be im- 
possible by the known general purpose algorithms. The incidence number of a face is 
the number of extreme rays contained in it; from our experience it is a good measure 
of the complexity of a face: in all polyhedral cones encountered by us so far, the 
facets with the highest incidence number are the ones of highest symmetry and they 
are also the ones whose dual description is the most difficult to compute. 

As a consequence, we begin the computation from the orbit with the lowest in- 
cidence number, since they are presumably easiest to treat and we may not need to 
treat all orbits, because of the following theorem due to Balinski: 

Theorem 3.1 ( jBa61f . see e.g. jZiQSi Theorem 3.14]) Let C he an m- dimensional, 
pointed polyhedral cone. Let G he the undirected graph whose vertices are the facets 
of C and whose edges are the ridges of C. Two vertices Ei,E2 are connected hy an 
edge F if Ei (1 E2 = F. Then, the graph G is {m — l)-connected, i.e., removal of any 
m — 2 vertices leaves it connected. 

Using the above theorem, we know that if the total number of facets in unfinished 
orbits is less than m — 1, then they cannot be adjacent to yet to be discovered facets 
and so we are done. In practice this simple criterion, which can be considered as an 
extension of Theoreme 7 and its corollaries in |Ja93] . is extremely useful and many 
difficult computations were finished by it. 

3.2 Recursive Adjacency Decomposition Method 

When there are only a few remaining facets to treat and we cannot apply Theo- 
rem [3]T1 our method is to use the Adjacency Decomposition Method on the remaining 
untreated facets of C recursively. The problem is that one may be confronted with 
a lot of cases to consider. In this section we present the method used to make this 
manageable in some cases. 

If a face F' satisfies F' < Fi <i C for a facet Fi, then there is exactly one other 
facet F2 such that F' <i F2 < C. Hence, if one applies the Adjacency Decomposition 
Method to Fi and F2, then one will compute the dual description of F' two times. 
The number of such repetitions increases as the recursion depth increases. 
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To handle this, we use a banking system. We store the representatives of orbits of 
{k — l)-dimensional faces (facets) F of a k-face F' with respect to a group Grp^/ C 
Auti?' of hnear automorphisms of F', which is not necessarily the full automorphism 
group of F' . 

If F' < F" and the Adjacency Decomposition Method is applied to F", then as a 
subtask one needs to find the orbits of facets of F' under the stabilizer Stab(Grpp^/,, F') 
of F' under the group Grp^,,. So, the problem is to obtain the list of facets under the 
action of Stab(Grpp,„, F') from a list of orbits under the action of Grp^,. 

When Stab(Grpp.//, F') is not a subgroup of Grp^./ (we did not require that Grp^./ = 
Auti?'), then we replace Grp^/ by the group generated by Grp^/ and Stab(Grp^//, F'). 
So, we can assume that Stab(Grp^//, F') C Grp^/. Take an orbit Grp^,/ F oi {k — 1)- 
dimensional faces of the fc-dimensional face F' . We find elements gi, . . . ,gr in Grp^,, 
such that 

r 

Grpp, = 1^ Stab(Grp^;/, F')5fj Stab(Grpp./, F), 

i=l 

i.e., we compute a decomposition of Grp^,, into double cosets (using a computer alge- 
bra system like GAP [GAPOSj l One then obtains Grp^,, F = IJ^^^ Stab(Grp^„, F^s-iF, 
i.e., the orbit Gipp, F splits into r orbits Stab(Grpp,», F')Fj with Fj = giF. This dou- 
ble coset decomposition is a classic enumeration technique, exposed for example in 
[B?00] and (Ke99] . 

To set-up a banking system as described, we need to be able to test if two poly- 
hedral cones are isomorphic and we need to compute their automorphism groups. 

3.3 Isomorphisms and automorphisms of polyhedral cones 

The symmetry group of a pointed polyhedral full dimensional cone in generated 
by extreme rays (ei)i<j<7v is the group of matrices A G GLm(M) such that there exists 
a permutation o" of {1, ... , A^} with Aci = erj{i) ioi 1 < i < N . This infinite group is 
called projective automorphism group (see |KaSc03j ). 

We are not aware of an algorithm to decide whether or not polyhedral cones 
are equivalent and to compute their projective automorphism groups. However, for 
the stronger notion of restricted equivalence introduced below, we can resolve those 
questions. Computing with a proper subgroup of the projective automorphism group 
is not a problem for the Recursive Adjacency Decomposition Method. However, it 
impacts the computing time since we may compute the dual description of projectively 
isomorphic polyhedral cones several times. 

Given a full-dimensional vector family (fj)i<i<Ar, we define the positive definite 
matrix 

N 
i=l 

Furthermore, denote by R the unique matrix R G S^q satisfying = R^ and 
{'Wi)i<i<N = {Rvi)i<i<N the image of the family (fj) under R. Define G{vi) to be the 
complete graph with vertices Vi and edge weights = ^ViQ~^Vj = {wi,Wj). 
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A restricted isomorphism of two vector families (fj)i<i<Ar and (fj')i<j<Ar is given 
by a matrix A such that there exists a permutation a satisfying Avi = v'^^--^ for 
t = l,...,N. 

Such a restricted isomorphism satisfies AQ^A = Q' with Q and Q' as above. One 
then checks that the matrix O = R'AR~^ is orthogonal and satisfies Owi = w'^^^y 
This implies Qj = c^^^^^^.^; i.e., the restricted isomorphism of the vector families {y.i) 
and (f •) corresponds to an isomorphism between the edge weighted graph G{yi) and 
G{v[). 

Now we need to prove that if a is an isomorphism between G{vi) and G(f^), then 
the equation Avi = f^^-^ admits only one solution. Clearly, we can assume a = Id 
and consider the equivalent equation Owi = w[. Since the Vi generate R™', we find a 
basis (vjj, . . . , Vi^) of M™. If P, respectively P', is the mxm matrix formed by (ifj^.), 
respectively {w'^J, then the equation Cij = c[j takes the form *PP = ^P'P'. So, the 
matrix O = P'P'^ is orthogonal and one has for any 1 < /c < m, 1 < j < A^: 

(w'i^^Owj) = {Owi^,Owj) = {wi^,Wj) = {wl,w'j). 

The above equation gives {w^, Owj — Wj) = 0; since the w'^ form a basis of M"^ we 
obtain the relation Owj = w'j, i.e., a restricted isomorphism between (vi) and (f •). In 
the same manner one proves that the restricted automorphism problem for a vector 
family (fi)i<i<iv is equivalent to the automorphism problem of the edge weighted 
graph G{vi). 

4 Implementation details 

The key part of our algorithm is to test if two {k — l)-dimensional faces of a k- 
dimensional face F are equivalent under a group Grp^;. of automorphisms of F. In 
practice, we represent {k — l)-dimensional faces by the set of indices of extreme rays 
contained in F and we use the command RepresentativeAction with the action OnSets 
of the GAP computer algebra system |GAP05j to test for equivalence. We compute 
the stabilizer of a face with the Stabilizer command. Building the full orbit, a strategy 
we could not consider, is used in [DFFSOl] and |DeDu03j . For some special groups and 
representations, like Sym(n) acting on n elements, it is very easy to test equivalence; 
this strategy is used in [An03j and |ChRo99j . 

In practice, the program nauty |MKa05j computes efficiently isomorphisms and 
automorphisms for vertex weighted graphs. Its "User's Guide" (Version 2.4, page 25) 
contains a description of a simple method to transform edge weighted problems into 
vertex weighted ones. 

5 Computational results 

To compute the dual description of the 10916 perfect domains we only had to use the 
Adjacency Decomposition Method in two cases: for Eg and for the Barnes-Coxeter 
lattice Ag (cf. |Ma03l Section 5.1]). It would also be necessary for the lattice Dg (cf. 
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[Ja92] ). but for D„ the list of neighboring domains is aheady known; see |RyBa79 



Theorem 15]. All other domains could be treated by the standard software Irs. 

Computing the dual description of Dom(QE8) (whose symmetry group is W{Eg)/'^^ 
of size 348364800) took several months of computer time. 

Most of the time of this computation was used for treating three orbits of facets 
of Dom((5Eg) generated by 66, 70 and 75 extreme rays respectively. The facet with 75 
extreme rays is interesting: its stabilizer under the action of the group W{Es)/'Z2 has 
size 23040, but when the automorphism group is computed, the size grows to 737280, 
therefore allowing us to finish the computation. We have no theoretical explanation 
for those additional symmetries. 

It is proved in |Wa71j that if a positive definite form A has || Min(A)| > 75, then 
A is arithmetically equivalent to a multiple of Qe^- The knowledge of facets of the 
perfect domain allows us to say a bit more: If || Min(A)| = 75, then A = Ai^i + A2A2 
with both Ai arithmetically equivalent to ^ind Aj > 0. If ||Min(A)| = 70, 
then A = XiAi + X2A2 with Aj > 0, Ai arithmetically equivalent to Qes ^2 
arithmetically equivalent to 

8 

It is also worthwhile to note that the form Qes is not contiguous to only two forms: 
Qas 8^nd the one with number 8190 (see http:/ /www. math. uni-magdeburg.de/lattice_geometry/| 
for the contiguities). 

The method developed here allowed us to treat the 8-dimensional case. Without 
using symmetry of perfect domains, one would be limited to dimension 6. Note 
that the algorithm developed by Jaquet in |Ja93j for dimension 7 is also a Recursive 
Adjacency Decomposition Method; he does not use a banking system and instead 
of using a standard software for the double description (which did not exist at that 
time) he uses a specific "Cascade algorithm" . We computed the facets of the perfect 
domain Dom(QE7) in less than a day. 

In dimension 9, a priori it seems currently very difficult to do a full computation 
of the list of perfect forms. Their number becomes even larger. We found more than 
500000 (see http://www.math.uni-magdeburg.de/lattice_geometry/) and there is no end 
in sight. Nevertheless, it is probably only a matter of time until the perfect forms in 
dimension 9 can be enumerated by Voronoi's algorithm. 
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6 Erratum to algorithm of section 2.3 



The algorithm given in section 2.3 is incorrect. It does not work when implemented 
and the computations of the paper were done with a different and correct algorithm, 
which is explained below: 

Input: Perfect form A and facet F of Dom(A). 
Output: Perfect form A' with Dom(yl') fl Dom(A) = F. 

(/,m)^(0,1) 

while A + uF ^ S^q or X{A + uF) = X{A) do 
if A + uF ^ 5^0 then 

u ^ {l + u)/2 
else 

(/, u) (m, 2u) 
end if 
end while 

while Mm{A + IF) C Min(g) do 

7^(/ + m)/2 

if X{A + 7F) > X{A) then 

else 

u ^ mm{{X{A) - A[v])/F[v] : v e Mm{A + -fF), F[v] < 0} U {7} 
end if 

if X{A + uF) = X{A) then 

I <— u 
end if 
end while 

The first section of the algorithm determines two values I and u which are lower and 
upper bounds for the value of the lifted parameter. The second section adjust until 
we find a correct value of the parameters based on shortest vector computations. The 
arrangement of tests is subtle but guarantees a termination in any case. 
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